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ABSTRACT 

We examine the radius evolution of close-in giant planets with a planet evolution model that couples the orbital-tidal and 
thermal evolution. For 45 transiting systems, we compute a large grid of cooling/contraction paths forward in time, starting 
from a large phase space of initial semi-major axes and eccentricities. Given observational constraints at the current time for a 
given planet (semi-major axis, eccentricity, and system age) we find possible evolutionary paths that match these constraints, and 
compare the calculated radii to observations. We find that tidal evolution has two effects. First, planets start their evolution at 
larger semi-major axis, allowing them to contract more efficiently at earlier times. Second, tidal heating can significantly inflate 
the radius when the orbit is being circularized, but this effect on the radius is short-lived thereafter. Often circularization of the 
orbit is proceeded by a long period while the semi-major axis slowly decreases. Some systems with previously unexplained large 
radii that we can reproduce with our coupled model are HAT-P-7, HAT-P-9, WASP- 10, and XO-4. This increases the number of 
planets for which we can match the radius from 24 (of 45) to as many as 35 for our standard case, but for some of these systems 
we are required to be viewing them at a special time around the era of current radius inflation. This is a concern for the viability of 
tidal inflation as a general mechanism to explain most inflated radii. Also, large initial eccentricities would have to be common. 
We also investigate the evolution of models that have a floor on the eccentricity, as may be due to a perturber. In this scenario 
we match the extremely large radius of WASP- 12b. This work may cast some doubt on our ability to accurately determine 
the interior heavy element enrichment of normal, non-inflated close-in planets, because of our dearth of knowledge about these 
planet's previous orbital-tidal histories. Finally, we find that the end state of most close-in planetary systems is disruption of the 
planet as it moves ever closer to its parent star. 

Subject headings: planetary systems - planets and satellites: general 



1. INTRODUCTION 

The precise mass and radius measurements for transiting 
exoplanets provide information about the planets' interior 
structure and composition, which are often apparently unlike 
that of Jupiter and Saturn. Indeed, it is the incredible di- 
versity of measured radii of transiting planets that has been 
most surprising. In the solar system, Jupiter and Saturn dif- 
fer in mass by a factor of three while their radii differ by 
only 18%. However, amongst exoplanets, planets with the 
same mass can differ in radius by a factor of two. A hope 
amongst planetary astrophysicists was that the measurement 
of the mass and radius, when compared to models, would 
cleanly yield information on planetary interior composition. 
Although there are clearly examples where this has been done 
successfull y, including heavy element rich plane ts such as HD 
149026b (iSato et al.ll200l iFortnev et al]l200fih and GJ 436b 
dGillon et al.ll2007^ 7in general modelers have been foiled by 
planets with very large radii, larger than can be accommo- 
dated by "standard" cooling/contraction models. 

Considerable work has been done in the past several years 
to understand the large radii of some planets, as well as the 
radius distribution of the planets as a whole. Explanations for 
the "anomalously" large planets have fallen into three cate- 
gories: those that are a current or recent additional internal 
energy s ource, which has stalled the interior cooling and con- 
traction (Bodenheimer et al. 2001; Guillot & Showman l2002l: 
Bodenheimer et all 12003b iGu et al.) 12003b IWinn & Holmanl 
2005 1 : iLiu et al.ll2008t Uackson et al.ll2008bt llbgui & Burrowsl 
2009), those that instead merely delay the contraction by 

1 Department of Astronomy and Astrophysics, University of California, 
Santa Cruz (neil@astro.ucsc.edu, jfortney@ucolick.org) 

- Lunar and Planetary Laboratory, University of Arizona 
(bj ackson @ lpl . arizona.edu) 



slowing the transport of i nterior energy (Bu rrows et alJl2007T: 
Chabrier & Baraffe 2007) , and those that invoke various evap- 
oratio n mechanisms (Baraffe et al. 200 41 lHansen & Barmanl 
2007). These are briefly reviewed in Fortneyl (120081) . 

Tidal heating as an ex planation for these l a rge-r adius plan- 
ets was suggested by iBodenheimer et all (1200 ll) for HD 
209458b a nd has been revisited frequently by other au- 
thors (e.g. IBodenheimer et al.l 120031 IWinn & Holmanl 120051: 
I Liu et al 112008b IGu et alJl2003M2004b1jackson et al.ll2008allbi: 
Ibgui & Burrows 200 91). We note that the me chanism of heat- 
ing by obliquity tides ( Winn & Holman 2005) has been cast in 
considerable double by several authors dLevrard et all 120071: 
iFabrvckv & Tremaine<2007l:|Peal3l2508b . 

At this time, tidal heating by orbit circularization is gener- 
ally believed to be the most important type. The largest un- 
certainties in the standard tidal theory is the "tidal Q" value, 
a standard parameterization of the rate of tidal effects. In 
this work, we use the standard notation for the planet tidal 
Q value as Q' p and the stellar tidal Q value as Q' s . Jupiter's 

< 2p value has been constrained to be between 10 5 and 10 6 
dGoldreich & SoterH l966). For tidal heating by circulariza- 
tion to take place, the planet must either initially have an ec- 
centric orbit or the system must be driving the eccentricity of 
the planet at recent times. 

The former scenario would have the following quali- 
tative stages. The planet is left w ith an eccentric or- 
bit through planet-planet interactions (|Rasio & Fordl Il996t 
IChatterjee et ani2008t iFord & Rasio1 l2008). Tides on the star 
gradually reduce the semi-major axis. These tidal effects ac- 
celerate as the semi-major axis decreases. Tides on the planet 
become more important and the planet's orbit circularizes; at 
the same time depositing orbital energy into the planet's in- 
terior. Scattering/tidal evolution models of this sort were re- 
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cently computed by Nagasawa et al. (2008). At this point, the 
system might be observed to h ave a fairly circul a r orb it and 
a larger-than-expected radius. Ilbgui & Burrows (2009|) use 
a coupled tidal-thermal evolution model, quite similar to the 
one we present here, to show that this scenario might be pos- 
sible for the HD 209458 system, and by extension, many hot 
Jupiter planets. Such a model is necessary to self-consistently 
explain a planet's radius in this picture. One potential issue 
with this scenario is that it can require large-radius planets to 
be observed at a "special time" since after the orbit is circu- 
larized, the planet may rapidly contract. 

Alternatively, some planets might be found in an equilib- 
rium state where their eccentricity is being forced by a third 
body while at the sa me time tides o n the planet are damp- 
ing the eccentricity (Mardling 2007). This is an attractive 
explanation because the planet might be found in an inflated 
state f or a long period of time. Previously, Bodenheimer et al. 
(120011) calculated the tidal power required to maintain the 
radius for HD 209458 b, Ups And b, and Tau Boo b, as a 
function of the assumed core size, in a stationary orbit. Re- 
cently, thermal evolution calculations with constant heating 
have been perfo rmed for TrES-4, XO-3b and HAT-P-lb by 
iLiu et alJ (120081) . who placed constraints on e 2 /Q' p - where e 
is the recent time-averaged eccentricity of the orbit. These 
calculations are useful for estimating the required recent tidal 
heating. In some cases, where the eccentricity is non-zero and 
a perturber is necessary to invoke, then this constant heating 
picture might accurately describe the recent thermal history 
of the planet. In many cases the eccentricity is observed to be 
close to zero, which either implies that a) the planet's eccen- 
tricity is at a non-zero equilibrium, but the planet's Q' p value 
is much smaller than inferred from Jupiter or b) the planet's 
orbit is circularized and this calculation does not apply. 

Clearly it is important to accurately measure the eccentric- 
ity of inflated systems to determine if either scenario is plau- 
sible. For many transiting systems, the eccentricity has been 
only weakly constrained with several radial velocity points 
and it is very difficult to distinguish a smal l eccentricity from 
one that is truly zero (Laughlin et al. 2005). For systems with 
an observed secondary eclipse, stronger upper limits on ec- 
centricity can be found based on the t iming of the eclipse 
dDeming et al.ll2Q05t ICharbonneau et al.ll2005l: iKnutson et all 
120091) . Note that secondary eclipse timing only constrains 
<?cos £1 so it is possible that some of these systems have much 
larger eccentricity, but it is unlikely. 

The above possibilities are also consistent with the popular 
planet formation and migration theories. These planets form 
while the protoplanetary disk is still present at much larger 
orbital distances an d migrate early in their life to small or- 
bital distances (e.g. lLin et al.lll996l) . After this initial phase, 
tidal evolution between the star and the planet occurs on Gyr 
time scales. The migration mechanism is important because 
it determines the initial orbital parameters for tidal evolution. 
There are multiple postulated migration mechanisms. 

1. Planet-disk interaction: Gravitational interactions be- 
tween the planet an d protoplaneta ry disk can exert 
torque on the planet (IWardl Il997alfbl) . These mecha- 
nisms tend to circularize the planet's orbit very early 
on and decrease the semi-major axis. The disk migra- 
tion time scales is signific antly shorter than the lif etime 
of the disk, as described in Papaloi zou et al.l (120071) and 
references therein. 



2. Planet-planet interaction: Gravitational interactions 
with other nearby planets can transfer orbital energy 
and angular momentum between the two bodies. This 
can result in quickly decreasing or increasing the orbital 
distance of one of the planets as well as producing non- 
zero i nitial eccentricity orbits. Using N-body simula- 
tions, dRasio & Ford 1996; Weidenschilling & Marzaril 
119961: IChatteriee et al.ll2008t iFord & Rasioll2008l) have 
shown that this effect can be important and can result 
in the inner bodies having initial eccentricity as large 
as 0.8, before tidal damping ensues. Other authors have 
investigated migration with coupled secular driving and 
tidal friction, which can operate on sim i lar timescales 
(IWu & Murray! 120031: iFaber et al.ll2005t IFord & Rasiol 
2006; Fabrycky & Tremaine 2007). There are also a 
handful of transiting planets that have non-zero eccen- 
tricity today, which can be explained by planet-planet 
interactions. It is also suggestive that many of these 
eccentric planets are more massive and have longer cir- 
cularization time scales. Since the circularization time 
is longer for massive planets, this observation is con- 
sistent with the idea that planets of all masses can have 
large initial eccentricity, but that the lower mass planets 
have circularized while the massive planets may still be 
circularizing. 

We expect that both of these mechanisms do happen to some 
extent. Therefore, we assume that a wide range of initial or- 
bital parameters are possible, and we following the orbital and 
structural evolution of planets from a wide range of possible 
initial eccentricities, as described below. In the absence of a 
theory to predict likely initial eccentricities for a given plan- 
etary system, we seek to understand the physics of the evolu- 
tion from a variety of initial states. 

Most of the detected transiting planets currently have small 
eccentricities consistent with zero. These can be explained by 
either migration mechanism. If the planet migrated through 
planet-disk interactions, then it would have zero eccentricity 
when tidal evolution began. If the planet migrated through 
planet-planet interactions, then the orbit may have circular- 
ized due to tides on the planet. 

2. MODEL: INTRODUCTION 

In this work we would like to test the possibility that tidal 
heating by orbit circularization can explain the transit radius 
observations for each particular system. A necessary condi- 
tion for this model is that a self-consistent evolution history 
can be found that agrees with all of the observed system pa- 
rameters. To check this condition, we forward-evolve a cou- 
pled tidal-thermal evolution model over a large grid of initial 
semi-major axis and eccentricity for each system. We perform 
this test for Q' s = 10 5 , Q' p = 10 5 and Q' s = 10 5 , Q' p = 10 6 5 . Also, 
for each system with non-zero current eccentricity, we emu- 
late an eccentricity driving source by performing runs with an 
eccentricity floor equal to the observed value. Later we also 
explore some higher Q' s cases. 

To properly understand the planet's thermal evolution, it is 
necessary to couple the planet thermal evolution to the orbital- 
tidal evolution. Generally planets with initial semi-major axis 
of 0.1 AU or less wi ll spiral into the star in Gyr timescales 
(Jackson et al. 2008a). This has a large impact on the incident 
flux on the planet, and therefore the loss of intrinsic lumi- 
nosity of the planet. For some systems, this more efficient 
cooling at early times makes it possible to achieve smaller 
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radii at the present. As the planet moves closer to the star, 
the tidal effects accelerate. If the orbit is eccentric, then at 
some point the planet's orbit undergoes a period of circular- 
ization. At this time a significant amount of orbital energy 
is deposited into the planet, which increases its radius. The 
question of this work is whether, at this stage, the system's ob- 
servables (a, e, age, R) can simultaneously be achieved in the 
model. After this stage, the p lanet may lose mass by Roche 
lobe overflow dGu et al.l20 03). which can temporarily prevent 
the planet from falling into the star. However, the planet's des- 
tiny i s to fall into the star dLevrard et al.ll2009t iJackson et al.l 
2009). These final stages of the planet's life, including the 
mass loss stage, are not modeled in this work. 

We typically find that tides on th e star are the dominan t 
source of semi-major axis evolution (IJackson et al.ll2008albh . 
When the eccentricity is large and damping, the tides on the 
planet can be the dominant semi-major axis damping source 
(IJackson et al.ll200"8at llbgui & Burrowsll2009l) . After survey- 
ing our suite of systems, we find that tidal heating can usually 
provide sufficient energy to inflate planetary radii as large as 
observed, but we do not always find an evolutionary history 
where the radius, semi-major axis, eccentricity and age all si- 
multaneously fall within the observed error bars. Regardless, 
we find that tidal processes are an important aspect of planet 
evolution, particularly for hot Jupiter systems. 

3. MODEL: IMPLEMENTATION 

The iFortney et al.1 (120071) gia nt planet thermal evo lution 
model has been coupled to the IJackson et al.l (l2008bl) tidal 
evolution model. Therefore, the semi-major axis, eccentric- 
ity, and radius of the planet all evolve simultaneously. The 
tidal power is assumed to be deposited uniformly into the en- 
velope of the planet. The planet structure model is assumed 
to be composed of four parts: 

1. a 50% rock/ 50% ice core (by mass ) with the ANEOS 
equations of state (Thompson 1990). The core does not 
participate in the the rmal evolution of the planet, as in 
IFortney et all (120071) . 

2. a H/He envelo pe with Y = 0-27, wh ich uses the equa- 
tion of state of Sau mon et al.l (119951) . The envelope is 
assumed to be fully convective and thus has constant 
specific entropy throughout. At each time step the en- 
velope is assumed to be in hydrostatic equilibrium. 

3. a series of radiative-convective, equilibrium chemistry, 
non-gr ey atmospher e mod els d escribed in more de - 
tail in IFortney et ail (120071) and IFortney et ail {2008). 
These grids are computed for the incident fluxes at 0.02, 
0.045, 0. 1 , and 1 AU from the Sun. This correctly deter- 
mines the atmospheric structure and luminosity of the 
planet as a function of the planet's surface gravity, in- 
cident flux from the host star, and interior specific en- 
tropy. In cases where the planet migrates to a semi- 
major axis with more incident flux than the innermost 
grid, then the boundary condition at the innermost grid 
is used. 

4. an extension of the atmosphere to a radius where the 
slant optical depth in a wide optical band (the Ke- 
pler bandpass) reaches unity. Therefore, all plotted 
radii are at t he "transit radius," as discussed by sev- 
eral authors (H ubbard et all 120011: iBaraffe et ail 120031: 



iBurrows et al.l l2003). The slant optical depth as a func- 
ti on of pressure is computed with the code desc ribed 
in iHubbard etaLl J200l and IFortney et al.1 d200ll) . We 
have found that the atmosphere height approximately 
follows the following relation 



h= 10 



8.74 1 eff 



(1) 



where h is the height in cm of the atmosphere from 
1 kbar (approximately the depth where the radia- 
tive/convective zone boundary lies) to 1 mbar (where 
the planet becomes optically thin), g is the planet's sur- 
face gravity (cgs), and T e g is the effective temperature 
in Kelvin. Taking into account this atmosphere height 
is significant when the planet has low grav ity or high 
effective temperature. In Fort nev et al.l (120071) . the radii 
at 1 bar were presented. 

The orbital -t idal evolution model is described in detail by 
IJackson et al.l d2008bl [2009) and references therein. The 
equations used in this work are 
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where a is the semi-major axis, e is the eccentricity, and P, 
is the tidal power deposited into the planet. This model at- 
tempts to describe tidal heating only by orbit circularization 
and ignores other forms of tides such as spin synchronization 
or obliquity tides, which are not believed to be as important. 
This model assumes that the star is rotating slowly relative 
to the orbit of the planet and is second order in eccentricity. 
Therefore the evolution histories that include periods when 
the orbit has high eccentricity should be regarded with cau- 
tion. Because there is a lot of other uncertainty with regard 
to tidal theory, we choose to us e this simple m odel instead of 
more complex models such as Wisdom (2008). For at least 
1 of the 45 systems, HAT-P-2, the planet-star system may 
be able to achieve a double tidally locked equilibrium state 
(star is tidally locked to the pl anet and the p l anet i s tidally 
locked to the star) as shown by Levrard et al. (2009); in this 
system it is not a good assumption that the s tar is rotating 
slower than the period of the orbit. However, iLevrard et alj 
(2009) find that this assumption is valid for most stars. We 
find that tidal heating is largest where e is not particularly 
large (< 0.4 falling towards zero) so this theory suffices for 
our purposes. Q' p is the tidal Q parameter of the planet and 
<2s is the tidal Q parameter of the star. In this work we have 
predominantly investigated cases when Q' p = Qi. = 10 5 as well 
as the case of Q' p = 10 6 5 , Q' s = 10 5 . Since th e Q value is in 
princ iple a function of the driving frequency dQgilvie & Linl 
2004), amplitude of the distortion, and internal structure of 
the body, the Q value for close-in extra solar giant planets is 
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potentially not equal to the Q value for Jupiter. If the Q value 
is a very "spiky" function of the driving frequency, then the 
system might spend a lot of time in a state where the tidal 
effects are occurring at a slow rate and quickly pass through 
states where tidal effects are rapid. The stellar Q value is typi- 
cally estimated through the observed circularization of binary 
stars orbits, but has also been estimated by m odeling the dis- 
sipation inside of a star (Ogi lvie & Linll2007h . 

We assume that the tidal power is uniformly deposited into 
the envelope of the planet. The net energy loss is given by the 
following equation: 

(L-P,)At = J TASdm. (5) 

where L is the luminosity at the planet's surface, At is some 
small nonzero time step, and S is the specific entropy. If 
P, > L, then the planet's envelope will be increasing in entropy 
and the planet's radius will increase. More typically, P, < L 
and the planet's entropy is decreasing and thus the planet is 
contracting. The power ratio P t /L is a useful measure of how 
important tidal effects are. It clearly indicates whether there 
is a net energy input (ratio larger than unity) or net energy loss 
(ratio smaller than unity). 

For a given radius, assumed core size and average inci- 
dent flux of the planet, R p oc -L„ et . Therefore, if we calculate 
Rnh, the radius contraction rate when there is no internal heat 
source, we can use the following relationship to calculate R p 
when there is an assumed P, tidal heating (or an input power 
of another source). 



Due to tidal migration the incident flux upon the planet in- 
creases with time. Based on the planet's incident flux at a 
given time, we interpolate in the 4 grids which include the in- 
cident flux level from the Sun at 0.02, 0.045, 0.1, and 1 AU. 
Here we neglect the more minor effect that parent star spectra 
can differ somewhat from that of the Sun. 

In order to examine all the plausible evolutionary tracks for 
each of the 45 transiting planets studied, we modeled their 
thermal evolution over a range of 

1. initial semi-major axis: the observed semi-major axis 
to five times the observed value. 

2. initial eccentricity: from to 0.8. 

3. core mass: 0, 10 M©, 30 M®, 100 M®. For very mas- 
sive planets we also consider core masses of 300 and 
1000 M ffi . Except for GJ 436b, HAT-P-llb, and HD 
149026b, the core was required to be at most 70% of 
the mass of the planet. For GJ 436b, we sample up to 
21 M ffi and for HAT-P-12, we sample up to 23 M e . 

Each of these possible evolution histories were run until ei- 
ther a) the time reached 14 Gyr, b) the entropy of the envelope 
became larger or smaller than the range of entropy values in 
the grid of hydrostatic equilibrium structures, or c) the planet 
reaches a small orbital distance ~ R s (realistically, the planet 
would be disrupted before this stage, but in this work we do 
not model the mass loss process). 

For each run, we searched the evolution history during the 
estimated system age range for times when the orbital param- 
eters were also within their observed range. If this occurred 



we then recorded the transit radius during these times and 
compared the range of achieved values to observed values. In 
situations where a good estimate on the age is not available, 
we searched within 1 to 5 Gyr. When a secondary eclipse 
constraint on the eccentricity is not available we assume that 
the eccentricity value is 0.025 ± 0.025 (i.e. the likely range 
is between and 0.05). In cases where the eccentricity is ob- 
served to be consistent with zero from a secondary eclipse, 
we assume that the eccentricity value is 0.005 ±0.005 (ie. the 
likely range is between and 0.01). We use the observed 
semi-major axis and error. We then search for instances of 
evolution histories during the possible age range that have an 
error-normalized distance less than 3 to the observed value. 
This distance is defined as 

yf (a, - a m ) 2 /a 2 a + (e,- - e m ) 2 /cr 2 (7) 

where a-, and e, are the orbital parameters for the instance 
of a particular run and a m , a a , e m , and a e are the mea- 
sured/assumed semi-major axis, semi-major axis sigma, ec- 
centricity, and eccentricity sigma. Planet orbital parameters, 
transit radii, and stellar parameters are from F. Pont's website 
at http://www.inscience.ch/transits/ and The Extrasolar Plan- 
ets Encyclopedia at http://exoplanet.eu/ 

4. GENERAL EXAMPLES 

Here we add different components of the model step-by- 
step, such that each effect can be appreciated independently. 
The two opposing effects of tidal evolution are late-time heat- 
ing that is associated with eccentricity damping and more ef- 
ficient early-time cooling due to initial semi-major axes that 
are larger then the present value. The four cases present are 
for a 1 Mj planet orbiting a 1 M Q star at 0.05 AU. In each of 
these cases we assume that the planet has a 10 M® core. 

Case 1 : no tidal effects, FigureQ] In the left panel, the solid 
line is the planet transit radius and the dot-dashed line is the 
radius at 1 kbar (near the convective-radiative boundary). In 
the right panel, the intrinsic planet luminosity is plotted as a 
function of time. As the planet contracts the luminosity of 
the planet significantly decreases. Without an internal heat 
source or semi-major axis evolution the planet's radius mono- 
tonically decreases with time. 

Case 2: no orbital evolution, constant interior heating, in 
Figure In this case the net output power is the difference 
between the intrinsic luminosity and a constant interior heat- 
ing source of unspecified origin. In these evolution runs, the 
planet stops contracting when the intrinsic luminosity is equal 
to the constant heating source. This is equivalent to when the 
ratio between the input power and the luminosity of the planet 
is equal to unity. The upper 3 evolution tracks (purple, cyan, 
and blue) all reach an equilibrium between the interior heating 
and luminosity of the planet within 2 Gyr, but the evolution 
runs with lower input power do not reach an equilibrium state 
in the 6 Gyr plotted. As expected, when there is more input 
power, the equilibrium radius is larger. In practice, the input 
power through tides or other processes will not be constant 
over gigayears, but a planet may be inflated to a radius such 
that it is in a temporary equilibrium state. 

Case 3: tidal orbital evolution, but without tidal heating, 
Figure [3] This case demonstrates how the orbital evolution 
due to tides effects the thermal evolution of the planet. Here 
we plot both the Q' p = 10 5 (tidal effects on the planet occur 
faster) and Q' p = 10 6 5 (tidal effects on the planet occur slower) 
cases with Q' s = 10 5 in black and red respectively. These 
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Fig. 1 . — : Radius and intrinsic planet luminosity evolution for 
a 1 Mj planet at 0.05 AU around a 1 M Q star without any tidal 
effects. In the left panel, the dashed line is the radius at 1 kbar, 
near the convective/radiative boundary at gigayear ages. The 
solid line is the radius where the atmosphere reaches 1 mbar - 
approximately the radius that would be observed in transit. 
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Fig. 2. — : Similar to FigureQ] but with various constant heat- 
ing applied in the interior of the planet. Moving from bottom 



to top, the constant heating rates are 10 , 10 , 10 , l(r ' and 
10 28 ergs" 1 . 

curves exactly track each other because the tides on the planet 
do not significantly contribute to the migration when the ec- 
centricity is small (here e = 0). When comparing Figure [JJto 
Figure [3] notice that in the second case, the power drops off 
more rapidly as the semi-major axis decreases. This is due 
to the increase in insolation by the parent star, which deepens 
the atmospheric radiat ive zone, lessening transport of energy 
from the interior (e.g. iGuillot et al.ll 1996b . Another result of 
moving the planet closer to the star is that there is an up-tick 
in the transit radius. This is due only to an increase in the 
effective temperature, which increases the atmosphere height. 
The semi-major axis evolution accelerates as the planet moves 
inward due to the tidal migration rate's strong dependence on 
semi-major axis. 

Case 4: tidal orbital evolution and tidal heating, Figure 
[4] We now put both the orbital evolution and corresponding 
tidal heating together. Black is the Q' p = 10 5 case and red is 
Q' = 10 6 ' 5 case. Notice that in the low Q 1 case, the planet cir- 




2 4 
me [Gyr] 



Fig. 3. — : Planet thermal evolution with orbit evolution, but 
without tidal heating. Transit radius, semi-major axis and 
the planet's intrinsic luminosity are plotted from left to right. 
Q' p = 10 s and Q' p = 10 6 5 cases are plotted in black and red 
respectively. 
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Fig. 4. — : Coupled planet thermal evolution and orbital evo- 
lution, gp = 10 5 and gp = 10 6 - 5 cases both with Q' s = 10 5 are 
plotted in black and red respectively. We plot, the radius evo- 
lution in the upper left, semi-major axis evolution in the upper 
right, ratio between tidal heating and intrinsic planet luminos- 
ity in the lower left, and eccentricity in the lower right. 

cularizes quickly and tidal heating becomes less important. In 
the high Q' p case, the planet is still undergoing circularization 
and significant tidal heating at late times. As a result, the ra- 
dius in the high Q' p case (slower rate of tidal effects in planet) 
can be larger than the low Q' p case (faster rate of tidal effects 
in planet) at late times. Both trials start out with fairly modest 
eccentricity (e = 0.3). 

In Figure [5] we compare the radius evolution in all four of 
these cases: Case 1 (no tidal effects, black), Case 2 (no orbital 
evolution, constant heating, blue), Case 3 (tidal orbital evolu- 
tion, but not tidal heating, red), and Case 4 (full tidal evolution 
model, cyan). The cases with tidal evolution are plotted for 
the high Q' p case. Clearly, when tidal heating is included (cyan 
or blue), it can result in a radius larger than achieved without 
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Fig. 5. — : Radius evolution in different cases. Cases 1, 2, 3, 4 
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including tidal heating (red or black). Since tidal heating is a 
time-varying quantity, the planet's radius when tidal heating 
will not be as simple as in Case 2. Generally, the planet will 
experience significant tidal heating when the orbit is being cir- 
cularized. At this time, the radius will increase, but after this 
time the radius of the planet will contract again. Also, because 
the planets in Case 3 (red) start at larger orbital distance than 
that of Case 1 (black), the radius contracts marginally faster 
when the planet is at larger semi-major axis. This is why the 
red line is lower than the black line before 2 Gyr. After this 
point, the transit radius increases in the red line case because 
the planet has moved close to the star, the effective temper- 
ature of the planet increases, and the atmosphere height also 
increases. 

To examine how different levels of internal heating affect 
the radius of the planet, we plot the planet radius after 5 Gyr 
as a function of mass in Figure [6] Again, these models as- 
sume a 10 core, at a orbital distance of 0.05 AU around 
a 1 Solar Mass star. In this figure, the black dotted line is the 
prediction of the thermal evolution model without tidal heat- 
ing. The red dashed line is the base of the atmosphere at 1 
kbar. Clearly, the height of the atmosphere is much larger for 
smaller planets due to their s maller gravities. The solid blue 
line is the radius relation from (Fort nev et al.ll2Q07h . The solid 
black lines are the radius of the planet given a constant heat- 
ing rate after 5 Gyr of evolution. The pink dotted curves are 
constructed in the same manner as the solid black curves, but 
required extrapolation (here, quadratic) off of the calculated 
atmosphere grid. At this point in time, most of these plan- 
ets have reached an equilibrium state where an equal amount 
of internal heating is balanced by the planet's intrinsic lumi- 
nosity. Clearly, the effect on the radius for a given heating is 
larger for smaller mass planets. 

5. RESULTS 

5.1. Specific Systems 

While we have computed the evolution history of 45 sys- 
tems, here we show representative calculations for particular 
samples of planets. These are TrES-lb, XO-4b, HD 209458b, 
and WASP- 12b, and are shown in Figures □ [8] [TOj andfTTI 
respectively. These four cases demonstrate qualitatively dif- 



Fig. 6. — : The dotted black line is the transit radius without 
any internal heating as a function of mass assuming a 10 M ffi 
core. In these models, we hold the planet at 0.05 AU around a 
1 Solar Mass star. The dashed red line is the 1 kbar radius — 
near the convecti ve/radiative zone bou ndary. The blue line is 
the relation from iFortnev et al] d2007l) . The solid black lines 
are the radius one would find if there were a constant heat- 
ing source (values between 10 24 and 10 29 erg s -1 ). The pink 
dotted lines were calculated in the same way, but required ex- 
trapolation (quadratic) off of the grid of atmosphere models. 

ferent cases. TrES-lb is a circularized planet with a "normal" 
radius value. XO-4b, HD 209458b, and WASP-12b are large- 
radii planets with a small relatively unconstrained eccentric- 
ity, zero eccentricity, and a nonzero value, respectively. In 
Figures [7] -QT| the transit radius evolution is plotted in the up- 
per left panel, the semi-major axis evolution is plotted in the 
upper right panel, the ratio between the tidal power and lumi- 
nosity is plotted in the lower left panel, and the eccentricity 
evolution is plotted in the lower right panel. The observed 
semi-major axis, eccentricity, and transit radius are plotted on 
each of the respective panels. The power ratio, tidal power to 
luminosity, describes how important tidal effects are to the en- 
ergy flow of the planet. When this ratio is somewhat smaller 
than unity, tidal heating is relatively un-important for the ther- 
mal evolution of the planet and when this ratio reaches or sur- 
passes unity, tidal heating plays a more significant role in the 
thermal evolution. In each of these figures, a set of runs were 
selected such that the orbital parameters and transit radius are 
closest to the observed values. 

TrES-lb is a transiting hot-jupiter planet with zero or small 
eccentricity and a typical radius observation. The system is 
composed of a 0.76 Mj planet orbiting a 0.89 M Q star with 
a 0.04 AU semi-major axis. Tidal heating is not necessary to 
invoke to explain this system; we demonstrate that this tidal 
model can still explain these kinds of modest radius systems. 
Possible evolution histories with tidal effects are shown in 
Figure [7] These possible histories are selected such that their 
orbital parameters at the current age agree with the observed 
values and the transit radius that is close to the observed value. 
We show various core sizes in different colors: black for zero 
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Fig. 7. — : Possible tidal/thermal evolution tracks for the 
planet around the star TrES-1. Black: no core. Red: 10 M ffi 
core. Blue: 30 core. Cyan dotted: 10 Me core evolution 
history without tidal effects. This is a 0.76 Mj planet orbiting 
a 0.89 Mo star. Upper left panel: transit radius evolution. Up- 
per right panel: semi-major axis evolution. Lower left panel: 
ratio between tidal power injected into the planet and intrinsic 
planet luminosity. Lower right panel: eccentricity evolution. 
Observed semi-major axis, eccentricity and observed radius 
are plotted in their respective panels. These evolution tracks 
were selected to have orbital parameters that agree with the 
observed values. Q' = 10 6 5 , Q' s = 10 s . 



core, red for a 10 M® core, and blue for a 30 M® core. The 
cyan dotted line is the evolution history of a non-tidal ther- 
mal evolution model with a 10 M© core. Notice the radius 
evolution of the non-tidal model doesn't differ significantly 
from the radius evolution of the corresponding 10 (red) 
tidal model. In these possible evolution histories with tidal ef- 
fects, the initial eccentricity is relatively small and tidal heat- 
ing doesn't dominate the energy flux budget (in the lower left 
panel, the power ratio is always less than 1). However, the 
orbit decays significantly due to tides raised on the star by the 
planet, which continues even at e = 0. These tides cause these 
planet to migrate from an initial semi-major axis of 0.05 AU 
to 0.04 AU with the assumed Q' s = 10 5 . Figure|7]demonstrates 
that this model easily explains the radius of TrES-lb with a 
core between 10 M® and 30 M® . 

There is a slight upturn in radius just before an age of 4 
Gyr. This is due to the heating of the planet's atmosphere at 
very small semi-major axis, and is not due to tidal power. As 
the planet reaches smaller orbital distances the incident flux 
it intercepts increases dramatically, leading to an enlarged at- 
mospheric extension, and greater tr ansit radius. This featur e 
is also present in the recent paper by Ibgui & Burr ows! (120091) . 
The tracks end when we stop following the evolution, with 
the assumption that the planet is disrupted or collides with 
the parent star. This is merely the first of many evolution 
tracks that we present with the end state being the disruption 
of the planet. This fi nding is essentially quite similar to that of 
iLevrard et al.l (|2009) who find that all of the known transiting 



planets, save HAT-P-2b, will eventually collide with their par- 
ent stars. Robust observ ational evidence for this mechanism 
was recently detailed bv lJackson et al.l (I2009I) . 

XO-4b is an inflated planet where the eccentricity has not 
been well constrained, du e to sparse radial velocity sampling 
dMcCullough et al] 120 08). In these cases we search for in- 
stances over the evolution histories where the eccentricity is 
between to 0.05, because we assume that a larger value 
would have been clearly noticed in radial velocity data. With 
this eccentricity constraint we show in Figure |8]that there is a 
narrow period of time when we can explain the inflated state 
with a recent circularization of the orbit that has deposited 
energy into the interior of the planet. The evolution curves 
shown here are for tidal parameters Q' p = 10 5 and Q' s = 10 s ; in 
the higher Q' p case, the radius evolution curves do not agree 
with the observed value. In Figure[8] we show black, red, and 
blue curves for evolution runs with no core, 10 M ffi core, and 
30 M ffi core respectively. The pink curve is an evolution his- 
tory for low initial eccentricity with a 30 M® core. Again, the 
cyan curve is a no-tidal evolution history with 10 M® core. 
Since tidal power is deposited mainly when the planet is be- 
ing circularized, high initial eccentricity orbits are required 
for these planets to experience significant later tidal inflation. 
Another interesting feature of this plot, is that when compar- 
ing the radius of the runs for different cores at any given time, 
we find that the radius is not always monotonically decreas- 
ing with core size. This shows that uncertain past orbital-tidal 
history can lead to uncertainly in derived structural parame- 
ters such as the core mass. 

As an example of the kind of calculation that was performed 
for every planet, in Figure [9] we show snapshots of the orbital 
parameters (a and e) of the ensemble of systems that are at 
some point consistent with the observed orbital parameters 
and age of XO-4b. Note that we do not require that the radius 
simultaneously also agree with the observed radius, but rather 
compare the range of possible radius values achieved by the 
model to the actual observed value. The black points are the 
original orbital parameters. The red points are the orbital pa- 
rameters for one of these runs at a later point in time (0.5 Gyr, 
1.5 Gyr, and 2.1 Gyr). The filled green circle marks the 1 a 
observed orbital parameters, while the dashed region is the 3 
a zone. 

HD 2094598b is a large-radius planet with ec centricity that 
has b een observed to be very close to zero dDeming et al.1 
120051) . The planet is observed to have a radius of 1.32 Rj 
and mass of 0.657 Mj. Therefore we require evolution histo- 
ries where the current eccentricity is < 0.01. Evolution his- 
tories for this system are shown in Figure [10] with Q' p = 10 5 

and Q' s = 10 5 . With these chosen Q values, we find that 
the planet could have experienced tidal heating at a previ- 
ous time, however by the time it has an eccentricity of 0.01 
or less the planet's radius has since deflated below the ob- 
served value. It is possible to find an evolution histories that 
agrees with t he observations by all owing different Q values, 
as shown by Ibgui & Burrows! (120091) . Although the tidal Q 
value is not strongly constrained and m ay even vary depend - 
ing on the configuration of the system dQgilvie & Linll2004b . 
it is our view that it makes the most sense to fix the Q value 
close to prior inferred values. Again, the black, red, and blue 
curves correspond to no core, 10 M® core, and 30 M® core 
sizes respectively. The cyan curve is a non-tidal thermal evo- 
lution history for a 10 M® core. In these cases, tidal power 
is sufficient to inflate the planet's radius to its observed value, 
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Fig. 8. — : Possible tidal/thermal evolution tracks for the 
planet around the star XO-4. This is a 1.72 Mj planet or- 
biting a 1.32 Mq star. Black: no core. Red: 10 Mq core. 
Blue: 30 M e core. Magenta: 30 M ffi core with a low initial 
eccentricity. Cyan dotted: 10 M ffi evolution history without 
tidal effects. Panels are analogous to Figure [7] The eccen- 
tricity that is marked in the lower right panel is our assumed 
possible range (0 to 0.05). These evolution tracks were se- 
lected to have orbital parameters that agree with the observed 
values. Q' p = Q[ = 10 5 . Notice that the tidal models initially 
have smaller radii than the non-tidal model because the tidal 
models are able to more efficiently cool at early times due to 
their larger semi-major axis. 

however we do not find evolution histories that also agrees 
with the other observed parameters — especially the eccentric- 
ity. In the semi-major axis evolution, there is a clear transition 
knee where the rate of orbital evolution decreases. The first 
phase is due to tidal effects of both the star and planet while 
the eccentricity is nonzero. The second phase is mainly due 
to tides on the star when the eccentricity is zero. 

WASP- 12b is a planet with an especially larg e radius of 
1.79 Rj with a non-zero eccentricity of 0.05 ( He bb et aD 
2009). An interesting property of this system is that t he planet 
is fill ing at least 80 % of its Roche lobe by radius dLi et al.l 
2009). Figure [TT1 shows evolution curves in black, red, and 
blue for no core, 10 M ffi core, and 30 M ffi core cases respec- 
tively when an eccentricity floor is imposed. Also, in cyan 
is the non-tidal model. In these tidal cases the tidal power 
increases in strength as the semi-major axis decays until the 
planet undergoes a rapid expansion. When the semi-major 
axis gets small enough, the tidal power exceeds the luminos- 
ity and the planet's radius rapidly increases. This happens 
both because the incident flux decreases the intrinsic lumi- 
nosity of the planet and tidal heating has a strong semi-major 
axis dependence (P t ~ a" 15 / 2 ). We do not model the mass 
loss process, w hich is likely to occur at late times for systems 
such as these dGu et al.l 12003) This should only be taken as 
evidence that if there was an eccent ricity driving com panion 
similar to mechanisms suggested by Mardling (2007), then it 
may be possible to heat this planet to quite large radii. 
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Fig. 9. — : Grid of evolution histories for XO-4b that were 
found to be consistent with the orbital parameters at a later 
time. These histories are not required to also have a radius 
value that is consistent with the observed value. These evo- 
lution runs assume a core size of 10 M ffi , Q' p = 10 5 , and Q' s 

= 10 5 . This serves as a sample for the type of calculation 
that was performed for every planet. Black: original orbital 
parameters of each run. Red: orbital parameters at a later 
marked time (0.5 Gyr, 1.5 Gyr, and 2.1 Gyr). The filled green 
circle is the 1 a zone, while the dashed region is the 3 a zone. 



5.2. Summary for Suite 

We have summarized our results for all 45 planetary sys- 
tems in Figures Q~2] and [13] for Q' p equal to 10 5 and 10 6 5 . 
In these figures, we have plotted the observed radius range 
(lower limit to upper limit) in black. The achieved radius 
range under various assumptions is plotted in color. Possi- 
ble radii are recorded in instances of the evolution histories 
when the orbital parameters and age all agree with the ob- 
served a, e, and age values (as defined previously, within 3 
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Fig. 10. — : Possible tidal/thermal evolution tracks for the 
planet around the star HD 209458. This is a 0.657 Mj planet 
orbiting a 1.101 M Q star. The planet has a radius of 1.32 Rj 
and an observed eccentricity of zero. Black: no core. Red: 
10 M ffi core. Blue: 30 M e core. Purple: 30 M e core with 
low initial eccentricity. Cyan dotted: 10 core evolution 
model without tidal effects. Panels are analogous to Figure|7] 

Q' P = Q' S = 10 5 . 
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Fig. 11.—: Possible tidal/thermal evolution for WASP-12b. 
This is a 1.41 Mj planet orbiting a 1.35 M Q star. The planet 
has a very large observed transit radius of 1.79 Rj and an ec- 
centricity of 0.05. In these evolution histories, we impose 
an eccentricity floor mimicking the effects of an eccentric- 
ity driving force. Black: no core. Red: 10M ffi core. Blue: 30 
Mq core. Cyan dotted: 10 core evolution history without 
tidal effects. Panels are analogous to Figure [7] Q' p = 10 5 and 

Gs = io 5 . 



error-normalized distance units of the observed value). The 
age of each system is often quite uncertain; since the possible 
radius values are sensitive to the age of the system, this is a 
large source of uncertainty for our results. For each planet, a 
range of radius values is plotted for up to five different suc- 
cessful types of models. These are models computed as dis- 
cussed in §|4] 

1. The full tidal evolution model is shown in purple. In 
this model the initial eccentricity was sampled from 
to 0.8 and the initial semi-major axis was sampled from 
the observed semi-major axis to 5 x the observed value. 
This is case 4 in §|4] 

2. The model with tidal migration but without heating is 
shown in green. We perform the same search procedure 
as in the full tidal model. This model is not meant to be 
physical, but to give us an understanding of how tidal 
orbital migration alone effects the planet's radius. This 
is case 3 in §|4] 

3. The "stationary" model is shown in blue with all 
tidal effects turned off. These are "standard" cool- 
ing/contraction mo dels, quite similar to those in 
iFortney et al.l (120071). These mode l s diffe r slightly than 
the models listed in IFortney et al.l (120071) in two ways. 
First, these models more accurately take into account 
the height of the atmosphere. Second, some of these 
models explore a wider range of core sizes. This is case 
lin&T] 

4. For planets whose current observed eccentricity is less 
than 0.4, the full tidal evolution with an maximum ini- 
tial eccentricity of 0.4 is plotted in orange. Because 
tidal heating in the planet is directly connected to ec- 
centricity damping, these runs serve as a demonstra- 
tion of relatively less tidal heating due to circulariza- 
tion. This is a subset of case 4 from Section |4] 

5. For systems where there is a measured non-zero eccen- 
tricity, we simulate the effects of an eccentricity source 
by performing the full tidal evolution with an eccen- 
tricity floor equal to the observed value. These cases 
are shown in red. This is essentially a combination of 
Case 4 and Case 2. 

For some planets, some of these "cases" were either not pos- 
sible to compute or in no instances were the observed param- 
eters consistent with the model parameters. For instance in 
cases when the observed eccentricity is larger than 0.4, the 
tidal evolution histories with 0.4 maximum initial eccentricity 
never are consistent with the observation. In these cases, no 
radius range is drawn. In some of the cases where tidal heat- 
ing is included, an evolution history is found where a large 
amount of energy is deposited into the planet while the orbital 
parameters are consistent with observations. These result in 
a maximum achieved radius that sometimes exceeds 2 Rj. In 
some of these cases, the planet will later cool off before the 
evolution stops. In other cases, the tidal power is sufficient to 
increase the planet's entropy beyond the maximum entropy of 
our grid, which ends the evolutionary calculation. In the fu- 
ture we plan to include mass loss and the subsequent evolution 
history. 

By comparing these models we find a few interesting pat- 
terns. When comparing the full tidal evolution model (purple) 
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Fig. 12. — : Observed planet radius (black) compared to a range of achieved model radii (colors) using £^, = 10 5 ;£^=10 5 . Planets 
are ordered by increasing incident flux according to their current observed parameters. Planets are marked with a * if they have 
nonzero observed eccentricity. The range of possible radius values under the full tidal evolution model is plotted in purple with 
initial eccentricity between and 0.8. The radius range for a model with tidal-orbital evolution, but without the tidal heating into 
the interior of the planet is plotted in green. The radius range for a standard stationary model without any tidal effects is plotted 
in blue. The radius range for the full tidal evolution model with a maximum initial eccentricity of 0.4 is plotted in orange. In 
cases where a nonzero eccentricity has been observed, the radius range with an eccentricity floor equal to the observed value is 
shown in red. 
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Fig. 13. — : Observed planet radius (black) compared to a range of viable model radii (colors) using Q' p = 10 6 5 ; Q' s = 10 5 . 
Qualitatively, we observe the same trends that were observed in Figure [T2al and ( 112bl i. A larger Q' p value decreases the rate of 
tidal effects via tides on the planet. Typically the tides on the planet from the star are responsible for circularizing the orbit, while 
tides on the star from the planet are responsible for decreasing the semi-major axis. In the larger Q' p case, the tidal circularization 
can be delayed for longer, which can make the possible radius of the planet larger. On the other hand, a larger Q' p also decreases 
the power deposited into the planet. 
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to the stationary model (blue), notice that there are some cases 
where the full tidal model has a larger maximum radius and 
other cases where the reverse is true. This can be understood 
to be caused by the two competing effects of tidal evolution. 
Tidal heating puts power into the planet and inflates the ra- 
dius, and tidal orbital evolution allows the planet to cool more 
efficiently at earlier times when the planet is less irradiated by 
the parent star. It is also useful to compare these two cases to 
the no heating model. The no heating model generally has a 
smaller maximum radius than the stationary model because of 
the second effect. The tidal model has a larger maximum ra- 
dius than the no heating model because of energy deposition 
into the planet. 

Often the model achieves large radius values through a 
recent circularization of an originally high eccentricity or- 
bit. During the circularization event (when the eccentricity 
drops significantly), tidal dissipation in the interior of the 
planet may deposit sufficient energy to significantly inflate the 
planet. The orange case (maximum initial eccentricity equal 
to 0.4) has been plotted to compare against the purple (initial 
eccentricity up to 0.8) to show how large initial eccentricity 
evolution histories contribute to the maximum achieved ra- 
dius. Note that in the low Q' p case in Figure Q~2] extremely 
large radii can be achieved for GJ 436b and HAT-P-1 1. This 
happens in our model through a recent rapid circularization of 
the orbit. 

It may also be possible to have tidal heating without large 
initial eccentricities if there is a eccentricity driving source in 
the system. In some cases, such as in WASP-6b or WASP- 
12b, the resulting tidal heating may be enough to explain the 
large transit radius. By comparing the red (tidal evolution 
with an eccentricity floor) to the purple (regular tidal evolu- 
tion), larger radius values can be achieved when the orbit is 
not allowed to circularize. 

Tidal evolution and heating clearly have important effects 
on a planet's evolution, but not all of the large-radius planets 
could be explained through this mechanism, given our chosen 
Q values. The planets HD 209458b, COROT-EXO-2b, HAT- 
P-9b, WASP- lb and TrES-4b have radii that are larger than 
achieved in our models in both the low and high Q' p cases. 
Typically, while it is possible to inflate the radius to the ob- 
served values, it difficult to find the system with an inflated 
radius and low current eccentricity. WASP- 12b was explained 
if we assume that its eccentricity is maintained. 

When comparing Figure [12] to Figure [13] it is interesting 
that some of the planets that are not explainable in the lower 
Q' p case can be explained with larger Q' p . Although Q' p = 10 5 
results in tidal heating being stronger than the Q' p = 10 6 5 case, 
it also results in circularization on a shorter time scale. In the 
Q' p = 10 6 ' 5 cases, it is often common for there to be a possible 
recent circularization of a high initial eccentricity orbit where 
no such history was found in the Q' p = 10 5 evolution runs. 

In Table 1, we have selected a set of the largest planets and 
listed various properties. In the left column, we list the ob- 
served parameters. For various core sizes, we list the achieved 
radius of the tidal model in the low Q' p and high Q' p cases, the 
estimated luminosity of the planet at its current radius, and the 
current contraction rate of the planet without internal heating 
(previously defined as Rnh)- Also, on the top row for each 
planet, we list the coefficient of tidal heating. This is defined 



as 




= — (GM.^M^fl- 15 / 2 x 10" 9 (9) 

This quantity allows one to get an order-of-magnitude idea of 
recent tidal heating given the more constrained properties of 
the system (radius of the planet, masses of the bodies, and 
semi-major axis). The actual tidal power will greatly depend 
on the eccentricity and Q values, which are more uncertain. 
The ratio between the luminosity of the planet and this coeffi- 
cient of tidal heating is a dimensionless number that describes 
how important tidal effects can be for a given core size. Cer- 
tainly, since P T oc e 2 and Q' p is quite uncertain, this ratio is 
not a strong test of tidal effects, but it is a simple way of test- 
ing how important tidal effects presently can be. Notice also 
that for an assumed tidal power, we can compute the present 
contraction rate using this table and Equation[6] 

When calculating the contraction rate, the planet is assumed 
to be located at the current observed semi-major axis, which 
determines the incident flux from the star, structure of the 
planet's atmosphere, and thus the intrinsic luminosity of the 
planet at each time. For these large-radius systems, the con- 
traction rate is often very fast. If we assume that tidal heating 
is the cause of large radii, but that an eccentricity driving com- 
panion is not present, then either the system is in a transient 
period or that this thermal evolution model is not correct. On 
the other hand, if we rule out transient explanations, then ei- 
ther a constant heating is present or it is necessary to invoke 
another mechanism. 

5.3. High Q' s cases 

Although Q[ is generally thought to be closer to 10 5 based 
on the observed circularization time in binaries, it is possible 
that that tidal dissipation in the stars is less efficient in the 
planet-star case. Since tidal evolution is not fully understood, 
the high Q' s case may or may not be physical. However, an 
advantage of this case is that it allows for orbital history solu- 
tions with a recent circularization. In this regime, the planet 
migrates inward at a slower rate and thus the circularization 
would occur at a later time. Also, after the tidal power is de- 
posited, the plan et is not rapidly m igrating into the star as in 
the low (2s cases. Ibgui & Burrows (2009) have suggested that 
high Q' s case can better explain the radius of HD 209458b. 

We have explored this parameter regime as shown in Table 

2 for five of the systems that we were not able to explain in the 
low Q' s cases. We test the cases Q' & = 10 6 and Q' s = 10 7 with 
both Q' p = 10 5 and Q' p = 10 6 5 . In the table the radius range is 
reported for a given core size, Q' p and Q[ model parameters, 
as well as the number of runs that were found at some point in 
time to be consistent with the observed age, semi-major axis 
and eccentricity of the system. 

Also, in Figure [14] we show snapshots in semi-major axis 
/ eccentricity space of possible evolution histories of HD 
209458 b that are consistent with the observed parameters. 
The black points are the original orbital parameters, while the 
red points are the orbital parameters at a later time. The green 
oval is the 1 a orbital parameters. The dashed green line is the 

3 <r orbital parameters, which we require an evolution histo- 
ries to fall within during the expected age range of the system. 
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Eccentricity was sampled from 0.2 to 0.8 in this particular 
case. 

We also show in Figure [impossible radius evolution histo- 
ries for the planets HD 209458b, WASP- lb, and CoRoT-Exo- 
2b. When Q' s is allowed to be larger, the qualitative effect is 
that the planet's semi-major axis decreases slower and thus 
the circularization event occurs at a later time. This makes it 
possible to sometimes achieve higher radius values at the ex- 
pected age of the system with the model. However, even for 
these high Q' s runs for these large-radius planets, only for two 
of the five can the observed radius be matched. 

6. DISCUSSION & CONCLUSIONS 

This paper presents a coupled tidal and thermal evolution 
model applicable to close-in extrasolar giant planets. The 
model is tested against 45 of the known transiting systems. 
Generally, tidal evolution yields two competing effects on the 
radii of close-in EGPs: 




0.03 0.04 0.05 0.06 0.07 0.08 
o [AU] 



(a) 



1 . Tidal evolution requires that, after planet formation and 
subsequent fast migration to a relatively close-in or- 
bit, the planet start at a larger semi-majo r axis than is 
currently observed (Jac kson et alj|200"8ah . This results 
in less incident flux at earlier times, which allows the 
planet to cool more efficiently and contract more at a 
young age, which moves the range of feasible model 
radii at the current time to smaller values. Generally 
this is a minor effect, but it is more important for cases 
when the current incident flux is larger. 

2. Tidal evolution deposits energy into the planet when the 
orbit is being circularized. This typically increases the 
radius of the planet at this time. If there is an eccentric- 
ity driving source for the inner planet, then tidal heating 
can be important for the duration of the planet's life. If 
the planet starts with a highly eccentric orbit, it might 
not circularize for gigayears. The semi-major axis of 
the planet's orbit will initially slowly decrease due to 
tides on the star. As the planet moves closer to the star, 
tides on the planet become more effective. This delay 
of circularization can sometimes allow tidal heating to 
significantly inflate planets multiple gigayears after for- 
mation despite these systems having shorter "circular- 
ization" time scales. 

We have shown that for the close-in giant planets that or- 
bital history can play a large role in determining the thermal 
evolution and current observed radius. While the effects are 
larger for planets with larger initially eccentricities, tidal evo- 
lution still affects the thermal evolution of planets with zero 
eccentricity as well. Varying amounts of time-dependent tidal 
heating are degenerate with the radius effects due to the core 
of a planet (or more generally, a heavy element enrichment). 

Since at the current time we are ignorant of the exact or- 
bital history, it is generally not possible to determine the mass 
of the core with complete confidence for any specific sys- 
tem. However, in cases when the radius of the planet is es- 
pecially small, a large core or increased heavy element abun- 
dance is required. For larger radius planets, it is not possible 
to determine the planet's core size because recent tidal heat- 
ing is degenerate with smaller core sizes. Furthermore, some 
systems likely have more complex orbital dynamics than de- 
scribed here due to the effects a third body. The uncertainty 
is increased since despite our expectation that tidal effects do 




0.03 0.04 0.05 0.06 0.07 0.08 
o [AU] 



(b) 




0.03 0.04 



0.07 0.08 



0.05 0.06 
a [AU] 

(c) 

Fig. 14. — : Grid of evolution histories (with initial e > 0.2) 
that were found to be consistent with the orbital parameters at 
a later time for the system HD 209458. These evolution runs 
assume there is no core, Q' p = 10 and Q' s = 10 6 . Black: orig- 
inal orbital parameters of each run. Red: orbital parameters at 
a later marked time (0.5 Gyr, 1 .5 Gyr, and 2. 1 Gyr). The filled 
green circle is the 1 a zone, while the dashed region is the 3 
a zone. 
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Transit Radius 

Q p =10 6 - 5 , Q,= 10 6 




o 



h- 

2 4 6 8 10 

time [Gyr] 

Fig. 15. — : Potential radius evolution histories for HD 
209458b, WASP- lb, and CoRoT-Exo-2b with no core, Q' p 

= 10 6 5 and Q' s = 10 6 (larger than our standard case). As usual, 
these evolution histories have been selected from an ensemble 
of possible initial conditions such that at some point during 
the estimated age of the system, the planet has orbital param- 
eters that are consistent with the observed values. 

occur, the rate that at which they occur (controlled by Q) is 
uncertain to an order of magnitude. 

This paper serves as a fo rward test of t he tidal theory for 
close-in EGPs outlined by I Jackson et alj (l2008bl) . who had 
previously only investigated heating rates backwards in time, 
from current small eccentricities from 0.001 to 0.03. Quite of- 
ten however, the forward modeling of these single-planet sys- 
tems, across a wide swath of initial a and e, is not consisten t 
with current eccentricities as large as Jackson et al. (2008b) 
assumed. If initial eccentricities were indeed large, then fi- 
nal circularization and tidal surge may indeed by fairly re- 
cent, but this cannot be expected to be the rule in these sys- 
tems. We have taken an agnostic view as to whether initial 
migration to within 0. 1 AU was via scattering or disk migra- 
tion. In_flie_Jormer, iiritial eccentricities up to 0.8 are possi- 
ble dChatteriee et al.l 120081) while in the latter the initial ec- 
centricity would be zero. The viability of tidal heating to ex- 
plain even some of the inflated planets with very small current 
eccentricities rests on the notion that planet scattering does 
occur, such that circularization (and radius inflation) can oc- 
cur at gigayear ages. The detection of misalignment between 
the planetary orbital plane axis and stellar rotation axis via 
the Rossiter-McLaughlin Effect (e.g. IWinn et~aT1l2007l 120081) 
is beg inning to shed light on migration. Fabrvcky & Winnl 
(2009) have found tentative evidence that is consistent with 
two modes of migration, one which may yield close alignment 
(perhaps from disk migration) and one with which may yield 
random alignment (perhaps from scattering), although to date 
only XO-3b in the pu blished literature shows a large misalign- 
ment (Hebrard et al. 2008). Further measurements will help 
to constrain the relative importance of these two modes of 
migration. 

Most of the systems investigated do not require tidal heating 
to match their radius, but these systems can also be readily 
explained when including tidal evolution. Some of the planets 
investigated can be matched with tidal heating that could not 
be explained with a standard contraction model. Depending 
on the Q' v value chosen, HAT-P-4, HAT-P-9, XO-4, HAT-P-6, 



OGLE-TR-211, WASP-4, WASP-12, TrES-3, HAT-P-7, and 
OGLE-TR-56 can all be explained with an evolution history 
with non-zero initial eccentricity. WASP-6 and WASP-12 can 
be explained by invoking a minimum eccentricity, which may 
suggest the presence of a companion. Other systems were 
not explained by the model for our chosen Q p values. This 
suggests that either Q' p and Q' s may be much different then 
our expectation or that other mechanisms are at work in these 
large-radius planets. 

This work should be taken as a simplified analysis of how 
tidal evolution can affect a planet's thermal evolution. Strong 
quantitative conclusions should not be drawn because of the 
large uncertainties in the tidal evolution model, especially at 
large eccentricity. Also, the rate of tidal effects may be a very 
strong function of frequency. If this is the case, the planet may 
spend a lot of time at certain states where tidal effects are slow 
and rapidly pass through states where tidal effects are more 
rapid. If a constant Q value can even be applied, the actual 
value is highly uncertain. The Q values that we choose were 
meant only to span the range that we considered to be likely. 
The rate of tidal effects may depend on the interior structure of 
the planet and may be different for different exoplanets. Also, 
this analysis only takes into account orbit-circularization tidal 
heating. 

The conclusion that should be drawn from this work is that 
a planet's tidal evolution history can play an important role 
on the planets' current radius, especially for systems that are 
born at semi-major axis less than 0.1 AU. In some cases, tidal 
heating could have inflated the radius of the planet in the re- 
cent past, even though tidal heating in the present might not 
be happening. In other cases, we were not able to explain the 
large-radius observations with our coupled tidal-thermal evo- 
lution model. This suggests that tidal heating will not be able 
to explain all of the l arge-radius planets, which has been a 
hope of some authors (Jackson et al. 2008b; Ibgui & Burrows! 
2009). For some of the planets that we are able to explain, 
we require a recent circularization, such that this model can 
only explain these observations if we at at a "special time" in 
its evolution. This has to be reconciled with the fraction of 
planets that have large radii that require such an explanation. 
Improved constraints on the eccentricities of these systems 
will better constrain recent tidal heating. 

A more robust treatment of the effects of tidal heating on 
transiting planet radius evolution may require a coupling of 
the model presented here to a scattering/disk migration model, 
which could derive the statistical likelihood of various initial 
orbital a and e configurations, which would then serve as the 
initial conditions to subsequent orbital-tidal and thermal evo- 
lution. This is important because for any particular planetary 
system the orbital evolutionary histor y of the close-in plane t 
may be difficult to ascertain. Recently Naga sawa et alJ (120081) 
have simulated the formation of hot Jupiters with a coupled 
scattering and tidal evolution code, and find a frequent oc- 
curence of hot Jupiter planets. A further coupled undertaking 
of this sort, to be compared with an statistically significant 
number of transiting planets, could be performed in the fu- 
ture. 
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TABLE 1: 

MODEL CALCULATIONS FOR SELECTED TRANSITING SYSTEMS 
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C T = 6.4 x 


10 23 




M p = 1.03 Mj 


0.0 


1.14- 1.38 


1.14- 1.22 


L = 5.0x 


10 26 


-6.7 x ur' 


R p = 1.36 Rj 


10.0 


1.12- 1.36 


1.12- 1.19 


L = 8.2x 


I0 26 


-i.i x ur 6 


a = 0.05 AU 


30.0 


1.08- 1.38 


1.07- 1.13 


L= 1.9 x 


I0 27 


-2.7 x IO" 6 


e = 0.00 


100.0 


0.93 - 1.10 


0.93 - 0.97 


L = 3.6 x 


8 


-5.8 x IO" 5 


WASP-1 








C T = 5.2 x 


r-*z 

10 26 




Mp = 0.87 My 


0.0 


1.16- 1.25 


1.16- 1.21 


L = 6.1 x 


10 2 " 


-1.4 x ur" 


R p = \MRj 


10.0 


1.13- 1.22 


1.13 - 1.18 


L= l.Ox 


I0 27 


-2.4 x ur 6 


a = 0.04 AU 


30.0 


1.07- 1.18 


1.07- 1.10 


L = 2.4x 


LO 27 


-6.2 x IO" 6 


e = 0.00 


100.0 


0.90-1.06 


0.90 - 0.92 


L = 5.1 x 


10 28 


-1.2 x IO" 4 


WASP-4 








C T = 1.3 x 


10 28 




Mp = 1.27 Mj 


0.0 


1.12- 1.20 


1.13-1.66 


L = 3.3x 


10 2 ' 


-3.8 x ur 6 


R p = \A$Rj 


10.0 


1.11 - 1.18 


1.10- 1.51 


L = 4.6x 


I0 27 


-5.4 x ur 6 


a = 0.02 AU 


30.0 


1.07- 1.11 


1.08-1.52 


L = 8.2x 


10 27 


-1.x 10~ 5 


e = 0.00 


100.0 


0.96-1.03 


0.96-1.18 


L = 7.5 x 


LO 28 


-8.7 x nr 5 


WASP-12 








Cr = l.lx 






Mp = 1.41 Mj 


0.0 




1.18-2.02 


L = 2.5 x 


io 28. 


-5.2 x ur 3 


Rp = U9 Rj 


10.0 




1.16-1.57 


L = 3.6x 


I0 28 


-7. x IO- 5 


a = 0.02 AU 


30.0 




1.12- 1.37 


L = 5.9x 


10 28 


-1.1 x IO" 4 


e = 0.05 


100.0 




1.01 - 1.11 









NOTE. — Various large-radius hot Jupiter planets have been listed. In the first column, we list the observed parameters of the system for 
reference. In the second column, we list an assumed core size. The achieved radius range for two different Q' p values is liested in the third 
and fourth columns. In the fifth column, we list relevant power quantities. The coefficient of tidal power is listed in the first row for each 
system. In the following rows, we list the luminosity of the planet for the assumed core mass. In the final row, we calculate Rnh, the radius 
derivative when there is no internal heating source. 
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TABLE 2: 

MODEL CALCULATIONS FOR SELECTED TRANSITING SYSTEMS 



System 


Core [Afe] 


Radius [Rj] (5,6) 


Radius [Rj] (5,7) 


Radius [Rj] (6.5,6) 


Radius [Rj] (6.5,7) 


HD209458 


M p = 0.69 Mj 
R P = \32Rj 


0.0 
10.0 


1.12- 
1.09- 


1.19(683) 
1.16(931 ) 


1.12- 
1.09- 


1.18 (737 ) 
1.15 (1136) 


1.15- 
1.12- 


1.32(816) 
1.27 (765 ) 


1.15- 
1.11 - 


1.31 (1036) 
1.25 (945 ) 


TrES-4 


M p = 0.93 Mj 
R p = \.l%Rj 


().() 
10.0 


1.16- 
1.13- 


1.22 (1291 ) 
1.19 (1285 ) 


1.16- 
1.13 - 


1.21 (849) 
1.18 (959) 


1.24- 
1.20- 


1.43 (665 ) 
1.37(512) 


1.19- 
1.16- 


1.37 (1205 ) 
1.33 (1154) 


HAT-P-8 


M p = 1.52 Mi 
R p = 1.58 Rj 


0.0 
10.0 


1.15- 
1.13- 


1.19 (1520) 
1.18(1515) 


1.15 - 
1.13- 


1.19 (1390) 
1.18 (1390) 


1.17 - 
1.16- 


1.28 ( 538 ) 
1.26(501 ) 


1.18- 
1.17- 


1.30(728) 
1.28(694) 


WASP-1 


M p = 0.87 Mi 
Rp = 1 .44 Rj 


0.0 
10.0 


1.17- 
1.14- 


1.21 ( 835 ) 
1.18 (829) 


1.17 
1.14- 


- 1.20(26) 
1.17 (297 ) 


1.23- 
1.20- 


1.48 (656) 
1.45 (636) 


1.19- 
1.16- 


1.39(1463 ) 
1.35 (1438 ) 


COROT-Exo-2 


Mp = 3.31 M/ 

Rp = \Al Rj 


0.0 
10.0 


1.12- 
1.11 - 


1.18 (1337 ) 
1.17 (1334) 


1.12- 
1.11 - 


1.19(1069) 
1.19 (1092) 


1.19- 
1.18- 


1.40(1243) 
1.39 (1242) 


1.13- 
1.12- 


1.33 (2127 ) 
1.32(2120) 



NOTE. — Achieved radius values for 5 systems with high Q' s for core size 0.0 and 10 M^ . The parameters used are denoted in the 
header with (log Q' p log Q' s ). In the body of the table, the range or achieved radius values is lested along with the number of runs found in 
parenthesis. 



